Nonideal quantum detectors in Bayesian formalism 
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The Bayesian formalism for a continuous measurement 
of solid-state qubits is derived for a model which takes into 
account several factors of the detector nonideality. In par- 
ticular, we consider additional classical output and backac- 
tion noises (with finite correlation), together with quantum- 
limited output and backaction noises, and take into account 
possible asymmetry of the detector coupling. The formalism 
is first derived for a single qubit and then generalized to the 
measurement of entangled qubits. 



I. INTRODUCTION 

The problem of continuous qubit measurement is 
of a significant importance for solid-state quantum 
computings because the measurement of a solid-state 
qubit typically requiEeSipsUsignificant time and thus can 
interplay nontriviallya'Bo&El with the intrinsic evolution 
of the qubit system. The evolution of a single solid-state 
qubit (without ensemble averaging) due to continuous 
measurement can be described by the Bayesian formal- 
ism (for review see RefE) which takes into account the 
noisy measurement output of the detector. The Bayesian 
formalism practically coincides with the version of the 
quantum trajectory formalism^ adapted to solid-state se- 
tups from the theory developed for quantum optics.tl 

One of the main predictions of the Bayesian formal- 
ism is the absence of the single qubit decohe*ence during 
the measurement by a good (ideal) detectorJlin contrast 
to decoherence of an ensemble of qubits. EM Moreover, 
the state of a solid-state qubit can be gradually puri- 
fied due to continuous measurement. In particular, this 
makes possible to monitor the phase of quantum coher- 
ent (Rabi) oscillations of the qubit. Such monitoring-cae, 
be naturally used in the quantum feedback controlE!rlL3 
of the Rabi oscillations which suppresses the qubit de- 
coherence due to environment (for q uant um feedback in 
quantum optics see, e.g., RefsJiallilolIa) . Another po- 
tentially useful application of the Bayesian formalism is 
a recent prediction that two qubits can be made fully en- 
tangled by their continuous measurement by an equally 
coupled detector .til 

The efficiency of the quantum feedback loop oper- 
ation crucially depends on the ideality (quantum effi- 
ciency) of the detector. For example, 100% synchroniza- 
tion between the qubit Rabi oscillations and desired puxa 
oscillations is possible only for 100% ideal detector.Ej 
Many other effects related to continuous measurement 
of solid-state qubits, which have been predicted using 



the Bayesian formalism (see, e.g., Refs.uOiiai 00) also 
depend significantly on the detector ideality. The ide- 
ality rj of a continuously operating solid-state detector 
can be generally defined as a ratio between the detector 
performance and the performance of a quantum-limited 
detector, in which the output and backaction noises are 
strictly related by the lower bound of an inequality sim- 
ilar to the Heisenberg uncertainty relation. More exact 
definition will be discussed later. 

A Quantum Point Contact (QPC) at low, temperature 
is theoretically an ideal (mantum detectoro that follows 
from the results of Refs.HH. A nearly ideal operaticm. 
of the QPC has been demonstrated experimentallyE3a 
The fact that a SQUID can theoretically reach the 
limit—of an ideal detector followstil from the results of 
Ref.Ej. A normal state single-electron transistor (SET) 
is not a good quantum detector at usual ooerating points 
above the Coulomb Blockade thresholdQO. However, 
its quantum .efficiency improves when we go closer to 
the thresholdEjEJ and becomes much better when the 
operating point is in the cotunncling range (below the 
threshold), in whinh-case the limit of an ideal detector 
can be achievedcU'Ea. Superconducting SET is gener- 
ally better than normal SET as a quantum-limited de- 
tector and can approach 100% ideality in the supercur- 
rent regimeo as welLas in the double Josephson-plus- 
quasiparticle regimeEl Finally, the resonant-tunneling 
setB can reach complete ideality in the small-bias lirait. 

In the simplest version of the Bayesian formalism!] a 
nonideal solid-state detectopis modeled as an ideal sym- 
metrically coupled detectorEil and a "pure dephasor" in 
parallel (environment or just extra backaction noise). In 
this case the nonideality leads to an extra term in the 
Bayesian equations, which introduces the gradual de- 
cay of the nondiagonal elements of the density matrix 
of the measured qubit. It was implied that such back- 
action dephasing is also equivalent to the extra noise at 
the detector output. However, the equivalence has never 
been proven explicitly, and this is one of the goals of the 
present paper. 

In a more advanced version of the Bayesian 
formalism,L3EJ a possible correlation between the out- 
put noise of a nonideal detector and the backaction noise 
is taken into account. However, the formulas for the evo- 
lution of the qubit density matrix in this case have been 
presented without any derivation, just from physical in- 
tuition. Moreover,, comparison of these formulas with 
the results of Ref .E2l for an ideal but asymmetrically cou- 
pled detector (which shifts the energy levels of the mea- 
sured qubit) reveals some difference. Even though the 
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difference is minor (second order in the detector response, 
which is assumed to be smaii), it points to some incor- 
rectness of the initial formulas of Ref.E3 (corrected formu- 
las can be found in RefE) . The main goal of this paper 
is to present a mathematical derivation of the Bayesian 
formalism for a nonideal detector with correlated out- 
put and backaction noises, using the phenomenological 
model which adds correlated classical noise to the quan- 
tum noise of an asymmetric ideal detector. We start 
with the measurement of one qubit and then generalize 
the formalism to the continuous measurement of an ar- 
bitrary number of entangled qubits. 

Notice that the issue of the asymmetric detector cou- 
pling to qubit has been recently discussed for a QPC 
in tefHjp-|©f, rtiie tunneling phase control by the qubit 
state.E3o'E3'E-3 For a small-transparency QPC the for- 
malism is significantly simplifiedcl] and is a direct gener- 
alization of the model of Reffl. We will use the results of 
Rcf.EHl to model an ideal asymmetric detector. 

While we model the detector nonideality by an addi- 
tional classical noise, let us-mention a different approach 
to the nonideality in Ref.EI in which a random fraction 
of electrons tunneled through the detector is assumed to 
be missing. In our opinion, such model is not well ap- 
plicable to solid-state detectors, even though it perfectly 
describes the inefficiency of a photodetector in a similar 
problem in quantum optics. 



II. MODEL 

We will use the phenomenological model of a nonideal 
solid-state detector of a qubit state shown in Fig. |l|. It 
consists of an ideal detector and three sources of addi- 
tional classical noise. We assume that the detector out- 
put is the noisy current I(t) (we have in mind a QPC or 
a SET as a detector) . The ideal detector is characterized 
by the output noise spectral density So [we assume flat 
("white") noise spectrum] and its backaction onto the 
measured qubit which will be called "quantum noise" . 
(Actually, because of the quantum relation between the 
output noise and the backaction, the output noise could 
also be called quantum; however, we will avoid such ter- 
minology, emphasizing the assumption that the quantum 
behavior does not propagate beyond the ideal detector.) 

The first source of an additional classical noise adds the 
noisy component with the white spectral density Si 
to the output Id(t) of the ideal detector, so that the final 
output is I(t) = Id(t) + £i(t). The second noise source 
is the classical noise £a(*) which is 100% correlated with 
(proportional to) the noise £i(i) and affects the qubit 
energy asymmetry e. [The qubit Hamiltonian is 

Kqb = | (4 C 2 - c\ci) + H (c\c 2 + cla), (1) 

where the tunneling strength H is assumed to be real 
without loss of generality] The relative magnitude of the 
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FIG. 1. Schematic of a nonideal solid state detector mea- 
suring a qubit state. The detector is modeled as an ideal 
quantum detector and three sources of additional classical 
noise: output noise £i(t) with the white spectral density Si, 
backaction noise &(t) fully correlated with and uncor- 

rected backaction noise ^(t). The total noise density So + Si 
of the output detector current I(t) includes the contribution 
So from the noise of ideal detector current Id(t). 



noise £2 (t) = A^i (t) is characterized by the parameter A. 
Finally, the third classical noise source is the white noise 
£3(i) which also affects the qubit energy asymmetry e [so 
that £ — > e + £2(1) + ^3 (i)] ■ The second and third noise 
sources together are obviously equivalent to one white 
noise source, partially correlated with (,i(t). However, 
we prefer to split it into the fully correlated and uncor- 
rected parts for clarity. Obviously, the qubit parameter 
H can also be affected by the detector noise; however, 
we do not take this effect into account, because the qubit 
dephasing is more naturally caused by the noise of its 
energy asymmetry e (which corresponds to the measured 
degree of freedom) and also because the induced noise 
of H is negligible, for example, for a single-Cooper-pair 
qubit measured by an SET. 

Let us start with the symmetric ideal detector, neglect 
all classical noises £1,2,3 if) and use the basic Bayesian 
formalism to describe the measurement process (i.e. the 
result of quantum backaction onto qubit) ;r-|rVpt the evo- 
lution of the qubit density matrix (t) isenJ Lil 

H 2 A/ 

P11 = -P22 = -2 — Impi2 + P11P22 —5— [h{t) ~ lo], (2) 
n Sq 

. s .H. 
P12 = 1 t P12 + 1 -r IP11 _ P22 ) 



-{p\\ - P22)^- [Id(t) - Jo] P12- 



(3) 



Here AI = I\ — I 2 is the detector response, I\ is the 
average detector current for the qubit state |1), I 2 is the 
average current for the state |2), and Jq = (I 1 +I»ll2 T , For 
the validity of Eqs. (|)-(§) we have to assumcu'Lrliil the 
weakly responding detector, |A/| <C Io, sufficiently large 
detector voltage (much larger than the qubit energies), 
and assume that the passage of individual electrons in the 
detector is much faster than the qubit evolution, Jo/e 3> 
{AH 2 + e 2 )/Ti, so that the current can be considered as 
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continuous. 

For simulations, Eqs. d)-@ should be complemented 
by the equation 



Id(t)-I = ^-(pU-p22) + to(t), 



(4) 



where £o (t) is the pure output noise of the ideal detector 
with flat spectral density So. 

Equations (@)~(@|) are written in the so-called 
Stratonovich formJiS which assumes symmetric defini- 
tion of the derivative, p(t) = lim r ^o[p(^ + r /2) — p(t — 
t/2)]/t. For the forward definition of derivative, £_(i)_ 
lim T _>o[p(i + t4i — P{t)]/ T (Ito form), the Eqs. 
transform intotil 



Pn = -P22 = -2 — Impi2 
n 



. e .H. 

Pl2 = 1 T Pl2 + 1 -r (Pll 

n, ft, 



2A7 „ , . 
P11P22 -5— ?o(t) , 



P22J 



,AJ (A/) 2 

"(Pll - ^22j-p-pi2 50(*J TE—P12 

Do 4D 



(5) 



(6) 



while the relation (Q) remains unchanged-, [The general 
rule of transformation is the following:E3 for an arbi- 
trary system of equations &i(t) — Gi(x,t) + Fi(x,t)£(t) 
in the Stratonovich form, the corresponding equations 
in the Ito form are ±i{t) = Gj(x, i) + i*i(x, t)£(t) + 
(S 5 /4) Y, k F k ( x > t) dF i ( x > *) / dx k ■} The advantage of the 
Ito form is that the averaging over the noise £o(t) is 
straightforward (we just need to eliminate terms with £ ), 
so it is easy to obtain the ensemble averaged evolution: 



Pn 



~P22 

e 



Pl2 = 1 T Pl2 

n 



-2 — Imp 12 , 
n 

(Pn - P22) 



(A/) 2 
4S 



■ P12 



(7) 
(8) 



On the other hand, the advantage of the Stratonovich 
form is the validity of usual calculus rules (which do not 
work in the Ito form) and therefore easier physical inter- 
pretation of equations. 

Let us emphasize that the single qubit in this model 
(which assumes ideal detector) does not decohere dur- 
ing the measurement process, as easier to see from Eqs. 
@ _ (||)- However, because of the probabilistic nature of 
quantum measurement, the ensemble of qubits does de- 
cohere (different qubits will go along different "trajec- 
tories"). The ensemble decoherence rate (AJ) 2 /4So is 
determined by this quantum randomness and therefore 
can be naturally called "quantum-limited" decoherence 
rate. (It can also be called "information-limited" ensem- 
ble decoherence, since its origin is the tendency of qubit 
state to evolve either into state |1) or |2), corresponding 
to the information acquired from the measurement O ) 

While it is not trivial to take into account additional 
classical noises £ 1 and £2 (this will be done in the follow- 
ing Sections), the account of the noise £3 is very simple. 
It leads to the additional dephasing term — 73P12 in Eqs. 



(||), (H), and (||), where 73 = S3 /Ah 2 is proportional to 
the spectral density S3 of £s(i). We will characterize this 
noise by the dephasing rate 73 instead of characterizing 
it by S3. (The relation 73 = S3 /4ft 2 can be easily derived 
adding £3 (t) to e in the Stratonovich form, then translat- 
ing the equation into Ito form, and averaging over £3.) 

The natural definition of the detector ideality factor rj 
in this case isQU'EJ 



V 



£0 

IV 



(9) 



where Tq = (A/) 2 /4So is the quantum- limited contribu- 
tion and Ts = To + 73 is the total ensemble dephasing 
rate. Simply speaking, this definition of ideality is the ra- 
tio between quantum contribution and total backaction 



III. IDEAL SYMMETRIC DETECTOR AND 
ADDITIONAL OUTPUT NOISE 

Let us now take into account additional output noise 
£i(t), while £2(4) is still zero. We also switch off £3^), 
since it is trivial to add its effect later. In order to derive 
Bayesian equations in this case, let us also assume H = 
e = ("frozen" qubit) and add the effects of H and 



later 



H 



0, Eqs. 



have a simple 



solutiorcrtLil which can be interpreted as a consequence of 
the "Quantum Bayes theorem" 



M T ) = 

P22{t) 
P12{t) 



P22W 



1 



1 

Mo) 



Pu( T ) 

M0)exp[-(ld- J 2 ) 2 t/Sq] 
Pii(0)exp[-(7J-/ 1 ) 2 r/So] 
Pii( T ) 

[pii(r)p 22 (r)] 1 / 2 



[Mo)Mo)F 2 

/ 3i 2 (0)exp[--^ j: ]exp[- 



(Id-Iofr -, 

s J 



pn(0)exp[ 



(/d-/l) 2 T l 

So -I 



P22(0)exp[ 



(ld-h) 2 T ] 

So J 



(10) 

(11) 

(12) 
(13) 



where Id is the average of the detector current during the 
time interval between and t: 



I d {t)dt. 



(14) 



5cw-|(|T^) is the consequence of the classical Bayes 
theoremt3 and Eq. 



Here Ec 

J_2|) says that the degree of the qubit 
purity is conserved. [It is easy to include the effect of 
finite e, which jus t leads to an extra factor exp(ier/?i) in 
Eqs. (y|) and fll3|); however, we will not do that in order 
to keep the formulas shorter.] 

Since the detector output is now I(t) = Id{t) + £i(t), 
we have to express Pij(r) in terms of I(t) and average it 
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over the noise £1 (t). Naively thinking, we have to use the 
substitution 



I d = I - x, I = - 



I(t)dt, x = - [ £x{t)dt, (15) 
o T Jo 



and average Eqs. ([!(]) and (|T^) over the noise 
contribution x using the weight factor p(x) = 
(2ttD 1 )- 1 ^ 2 exp(-x 2 /2Di) where D x = Si/2t is the vari- 
ance of x. However, this is not a correct procedure be- 
cause the probability distribution of x is correlated with 
I (though it is not correlated with Id). So instead, we 
have to use the conditional distribution of x for a given 
7: 



(16) 



p(x) = P{x)l J P{x')dx' 
cxp(-a; 2 /2L>i) 



P(x) = 



Pii(O) 



exp[-(I-x-h) 2 /2D ] 



+/O 22 (0) 



exp[-(I-x-l2) 2 /2D ] 



V2TD~ Q 



(17) 



where D = Sq/2t. (Let us stress again that both Id and 
x are assumed to be classical quantities, so the derivation 
is not applicable to the case when the detector output Id 
is involved in further quantum interactions.) Introducing 
the weight factor p(x) into Eq. ([To|), substituting I d = I— 
x, and integrating over x, we get the averaged equation 



/fn(r) 



p 22 (0) exp[-(7- J 2 ) 2 t/S s ] 
Pii(O) exp[-(7-/ 1 ) 2 r/5 s ] 



(18) 



where = Sq + Si is the total output noise. The only 
difference compared with Eq. ([lO]) is the change of Id into 
/ and change of Sq into Ss (this is quite expected since 
pn behaves as a classical probability and the classical 
Bayes formula still works). 

To calculate Pi 2 (t) averaged over the noise we 
have to do a similar procedure. We multiple Eq. (13) by 
the weight factor p(x), use substitution Id = I — x, and 
integrate over x. In this way we obtain 



pi 2 (0)exp[ 



(AJ) 2 - 
4S 



|exp[ 



pii(0)exp[ 



(I-h) 2 T ] 

s E -I 



p 22 (0)exp[ 



(19) 



Comparing Eqs. (|l^ ) and (|l9|), we see that Id changes 
into / and So changes into Sh, except in the second fac- 
tor of the numerator, where So remains unchanged. Let 
us represent this factor as exp[— (A/) 2 T/4Sy exp(— Jit), 
where 



7i 



(AI) 2 S 1 



(20) 



So, the effect of additional output noise £i on the Eqs. 
(|l0|)-([r3|) is the following: the output current Id from 



the ideal part of the detector changes into the output 
current /, the spectral density So corresponding to the 
ideal part of the detector changes into the total output 
noise Se, and the nondiagonal matrix element acquires 
the dephasing factor exp(— 7it). Differentiating the new 
equations over time (if we do it in a simple first-order 
way, we automatically get equations in the Stratonovich 
form) and adding terms due to H, e, and noise £3, we 
obtain 



-R22 

e 



Pl2 = 1 T Pl2 

n 



-2 — Im p 12 

n 

AI 



2AI 

•pllp22-5-[/(t)--fo], (21) 



P22) 



-{pit - P22)-^~ [I(t) ~ h\ P12 
<3E 



Ss 



(71 + 73)P 



12- 



(22) 



We see that the effect of the extra output noise S,i{t) is 
similar to the effect of the extra backaction noise ^(t) 
and both lead to the qubit dephasing. 

From physical reasoning, the way of separation of the 
detector into the ideal part and additional noise sources 
£i(i) and £3^) is arbitrary, as long as the total output 
noise Ss and total ensemble qubit dephasing rate are 
fixed (in other words, Ss and are the only physically 
relevant quantities). It is easy to check that Eqs. (I2T 
(|22j) satisfy this requirement because 



71 + 73 = r 5 



(All 
4S^ 



(Alf 
4S 



- 73- 



(23) 



[The total ensemble dephasing rate Ts can be formally 
found from Eqs. (21 )-(|22j) by translating them into Ito 



form that adds ensemble dephasing rate (A/) 2 /4Ss).] 

Comparing Eqs. @-(||) with Eqs. we natu- 

rally introduce a more general definition of the detector 
ideality: 



V 



(A/) 2 /4S s 



(24) 



which is again the ratio of the quantum-limited part of 
the ensemble qubit dephasing and its total dephasing 
rate. (Notice that the numerator is not the "real" quan- 
tum backaction determined by So, but the information- 
limited backaction determined by Ss-) In particular, for 
our model in the case £3 (t) = (no classical backaction) 
we obtain 77 = So/ (So + Si). 



IV. CORRELATED OUTPUT AND 
BACKACTION NOISES 

Now let us add the classical backaction noise £ 2 (t) 
which affects the qubit energy e [so that e — > e + £ 2 (i)], 
and which is 100% correlated with the output noise 
source, £ 2 (£) = A£i(t). We again start with Eqs. (|h])- 
( |l3| ) for the "frozen" qubit. Addition of £ 2 (t) does not 
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affect diagonal matrix elements, so Eq. (|T^) remains un- 
changed. In order to calculate Pi2(t) averaged over the 
noises £1 and £2, we multiply Eq. (|l3|) by the factor 
exp[i?i _1 ^2(t)dt] = exp(iAxT/H) and average it over 
x with the weight p(x) given by Eqs. (]T^)~(|l7|). This 
leads to the equation 



pMt) 



pi 2 (0)exp[ 



pn(0)exp[ 



(i-h) 2 



P22(0)exp(- 



(i-h) 



xexpI-U-IoM-rlexph^-^r] 



(25) 



in which only the second line is different from Eq. (|19|). 

Differentiating this equation over time (again, we au- 
tomatically obtain the Stratonovich form) and adding 
terms due to H, e and £3, we get 



where 



(A/) 2 
45 s 



K 2 St 



(32) 



Eqs. (|30|)-(|32|) is the main result of this Section. Com- 
paring them with-jSisnilar equations presented (but not 
derived) in Refs.E30, we notice a difference: the .term 
iK[I(t) — Io]/9i2 was incorrectly replaced in Refs.EdO by 
the term iK[I(t) — {pnh + P22h)]pi2- Notice though 
that the effect of their difference iK(AI /2)(pn — ^22)^12 
is minor since pn — P22 is the oscillating magnitude and 
averages to zero. As will be mentioned later, Eqs. j30|)- 
( |3~l| ) in the case 7 = (this is possible for asymmetric 
ideal detector) coincide with the corresponding equations 
of Ref.E3 if e includes the detector-induced shift. 

To translate Eqs. ©-© fr om Stratonovich to Ito 
form, notice that 



. £ . H , 

Pl2 = 1 T Pl2 + 1 -T" (Pll 

n n 
AI 



i S 

P22) + ^[I(t)-I ]A^Pl2 



"(pil - P22)-£T [I(t) - I ] P12 ~ (71 + 72 + 73)Pl2, (26) 



where 71 is given by Eq. ( p0[ ) and 72 is 



72 



SqS\ A 

"SrTJh 2 



(27) 



Physically relevant parameters of the detector are the 
total output noise 5s = So+S±, total ensemble dephasing 
rate r E , and the correlation between output and backac- 
tion noises. Since three sources of the backaction noise 
in Fig. [l] are uncorrelated, the ensemble dephasing rate 
is 



/(*) - ^0 



AI 



(pll -p22)+6> (*)+&(*) 



(33) 



and the sum of two output white noises t;o+i(t) = £o(t) + 
£i(t) is the white noise with the spectral density 5s- 
Then using the standard rule of translation, we obtain 
Ito equations 



P11 

Pl2 



-P22 



i Jpi2 



H 

-2 — Impi2 

h 



AI 



2AI , . . 
P11P22 fo+i W, 

P22) +iK£ 0+ i(t) P12 



(34) 



(Pll - ^2)^7 &+1 (t) pi2 

. (AI) 2 K 2 Sx\ 

7 + W + ^J P12 ' 



(35) 



T s = (A/) 2 /45 + A 2 S 1 /4h 2 



73 



(28) 



(the same result can be obtained by translating Eq. ( |2q ) 
into Ito forip_£LLid performing ensemble averaging). Fol- 
lowing Refs.E3'EJ, let us characterize the noise correlation 
by the magnitude (real number) 



K 



S E i 



(29) 



where S £ i is the mutual spectral density of the detec- 
tor output noise and induced fluctuations of e (strictly 
speaking, S e j in our notation is only the real part of the 
mutual spectral density, while the imaginary part can 
formally describe the detector responsecj). In our case 
K = AS\/KS-£. Expressing qubit evolution only in terms 
of physically relevant quantities, we obtain 

H 2AI . , . . . s 

P11 = -P22 = -2 — Impi2 + P11P22 -5— [I(t) - ioj, (30) 
a 5s 

s H 

P12 = i t P12 + i -z- (P11 - P22) + iK[I{t) - I Q ] P12 
n n 

-(Pll - P22)^" [I(t) - Io] P12 - 7P12, (31) 



while the relation between output current I(t) and pure 
noise £o+i (t) is still given by Eq. (|33|). (Notice that 
the above mentioned incorrect term is correct in the Ito 
form of the equation, so the mistake was due to mixing 
of the Stratonovich and Ito forms.) The corresponding 
ensemble-averaged equations can be obtained by erasing 
terms containing £o+i(*) m Eqs. (|3^)-(|35|): 



P11 = -P22 = -2 — Impi2 , 
n 



P12 



i|pi 2 



i -r (pll - P22) - T E pi2 



(36) 
(37) 



where the total ensemble decoherence rate is given by 
Eqs. (H) and/or (§§). 

Since 7 > [otherwise solution of Eqs. (|3^)-(|35|) would 
violate inequality |pi2| 2 < P11P22], we obtain the funda- 
mental limitation for the ensemble dephasing: 



Ty, > 



(All 
45 s 



K 2 Sy 



(38) 



Besides the definition of the detector ideality 77 given by 
Eq. (||) which would give 77 = 1- (7 + if 2 5s/4)/T E , it-is 
natural to introduce another definition of the idealityO 
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J_ 



(AT) 2 /4S S + K 2 S^/A 



(39) 



since the term K 2 S^/A does not correspond to the de- 
phasing of a single qubit. One more possible, definition 
of ideality (which also gives 100% if 7 = 0) isQ 

(AT) 2 /4S S 



fj 2 = t = ^ (40) 

12 l + 7 /[(AT) 2 /4S s ] r s -X 2 S s /4' ^ 

so that (^2) ^ 2 directly correspondsS to the total energy 
sensitivity of the detector in units of h/2. In the case 
K = all the definitions coincide, 77 = fj = 



V. ACCOUNT OF ASYMMETRIC IDEAL 
DETECTOR 

So far we have assumed that the ideal part of the detec- 
tor in Fig. [ijdoes not induce the shift of the qubit energy 
asymmetry e (i.e. in our terminology assumed symmet- 
rically coupled detector). However, in general the cou- 
pling with detector changes e, so it should be treated 
self-consistently.y As an example, the operating point 
of an SET as a detector is slightly shifted by different 
charge states of a measured qubit. This generally leads 
to the change of the average potential v of the SET is- 
land, which affects back the qubit energy asymmetry e. 
Notice that v can also be temporary shifted by a fluc- 
tuation of the current through SET, leading to the cor- 
relation between the output and backaction noises. So, 
in this example the shift of e and noise correlation are 
closely related. Similar situation occurs when as a detec- 
tor we use a QPC, which legation relative to the qubit is 
geometrically asymmetric. ffil Then the qubit state af- 
fects the phase of the QPC current (of course, this should 
happen before the current becomes a classical quantity) , 
and in return each electron passing through the QPC af- 
fects the phase difference between qubit states |1) and 
|2), thus leading to effective shift of e. Correspondingly, 
the noise of the QPC current causes the correlated noise 
of e, so again these effects are closely related. 

The asymmetric coupling can be relatively easy taken 
into account for a smalL-transparency QPC using the 
model analyzed in RefHl The detector and its inter- 
action with the qubit are described by Hamiltonians 

Ti-det = Eiajai + E r a\a r + ^(Ta^a/ + T*aja r ), 

l r l } r 

(AT 1 AT 1 * \ 
-— aJaH —aja r j, (41) 

in which we neglect the dependence of tunnel matrix el- 
ements T and AT on the electron states (l,r) in elec- 
trodes. n The only difference of this model from the model 
of Ref.0 [which leads to Eqs. is the possibility 

of complex T and AT (actually, T can be assumed real 



without loss o£-generality) . Following .the procedure de- 
veloped in Ref.0, it is possible to showES that in the cor- 
responding Bloch equation for p™ 2 (where is the den- 
sity matrix with account of the number n of electrons 
passed through the detector) the termn y/ Iihe -1 Pi? 1 
should be replaced by exp(i(/j)v / ^i^2e~ 1 p™ 2 _1 ' where ip = 
arg[(T + AT/2)(T* - AT*/2)]. Therefore, each electron 
tunneling through the detector adds the phase ip to pi2, 
thus affecting the qubit energy asymmetry e. 

The assumption of weak detector response implies 
I AT I < |T|, so that \ip\ < 1.E3 The extra phase leads 
to the extra term i(9/h)IdPi2 in Eq. (|J) where 9 = Tup/e. 
Separating it into the average and fluctuating parts, we 
obtain the following equations for the asymmetric ideal 
detector in the Stratonovich form: 

TT 2 AT 

P11 = -P22 = -2 — Impi2 + P11P22 -5— [Ttf(i) - T ], (42) 
h bo 

s TT 9 

P12 = i t P12 + i -t" (pn - P22) + i t [hit) - T ] P12 
n n n 



-{pn - P22)^~ [Id(t) - T ] pi2, 

>~>0 



(43) 



where e = e + Ae and Ae = 9Iq = hiplo/e. Equations 
(^2])-(^3|) in ltd form have beeruqbtained in Ref.o (no- 
tice a different definition of 9)tB Let us mention that 
even though Eqs. (42)-j4^) have been derived for a par- 
ticular model of the detector [Eq. (|4j])], it is expected 
that they are applicable to a significantly broader class 
of asymmetric ideal detectors (i.e. finite-transparency 
QPCs, quantum-limited dc SQUIDs, etc.). Notice that 
Eqs. ©-(p|) are formally similar to Eqs. ©-@ ex- 
cept 7 = (ideal detector) and s is replaced by the self- 
consistent value e. 

The solution of Eqs. (ph-(p3|) in the simple case H = 
is still given by Eqs. (|10[)-(|ll|Jfor pn and P22, while Eqs. 
(|l2^-(|l3|) for pi2 should be replaced by 

^ = ^ TtS ^ e^~^ (44) 



/ n x r (AI) 2 T (Id-Io) 2 Tl 



So 



exp[i gllgCSzjojr ] 



pn(0) exp[ 



(Jd-/l) 2 r i 
So J 



P22(0)exp[ 



(Id-l2) 2 T ] 
So J 



(45) 



Let us now add the classical noises £1, £ 2 and £3 (see 
Fig. |l|) to our model of asymmetric ideal detector. Us- 
ing the procedure explained in two previous Sections, 
we multiply Eq. (45) by the factor exp(L4xT/7i) where 
x = t -1 f Q £i(t)dt, average the resulting equation and 
Eq. ( [l0| ) over the distribution p(x) given by Eqs. (|l6|)- 
( |l7| ) , differentiate the result over time, and add the terms 
due to H and £3. In this way we obtain the following 
Stratonovich equations: 

H 2 AT 

Pn = -P22 = -2 — Impi2 + P11P22 s— [^(*) _ to], (46) 
Ti Sy, 

. i . H , 
Pu = 1 -g P12 + 1 -j- (Pn - P22 ) 



G 



i ASi + 9S 



[I(t) - Io] P12 



-{pn ~ P22)^~ [I(t) - h] Pu 



where 



7 



7P12, 



- 73- 



(47) 



(48) 



4S S S S E 4?i 2 
It is interesting to notice that the single-qubit decoher- 
ence rate 7 can be decreased by adding the backaction 
noise £2^) = ^Ci(^) with A = 9 (i.e. having the same 
correlation with the output noise as the quantum back- 
action) , even though it increases the ensemble dephasing 
rate 

r s = (Ai) 2 /4S + 9 2 s /^n 2 + A 2 s 1 /m 2 + 73 . (49) 

Introducing the total correlation between the output 
noise and the backaction (including quantum) noise, 

ASi + 9S 



K = 



(50) 



we reduce Eqs. ©-© to Eqs. ©-© with the only 
substitution of e by the self-consistent value e. There- 
fore, the corresponding Ito equations are still given by 
Eqs. (|34|)~(|35|) anc j the ensemble-averaged equations are 
still given by Eqs. (|36|)-(37) (with e — > e), while the 



relation between the single-qubit decoherence 7 and en- 
semble decoherence still satisfies Eq. (|3^). Thus, we 
conclude that the Bayesian description of the measure- 
ment process given by Eqs. (|o|)-((3^), which is expressed 
in terms of the physically relevant quantities K, Ss, and 
Ts, remains valid in the case when the ideal part of the 
detector is assumed asymmetrically coupled to the qubit. 
Similarly, the limitation (38) for remains valid and the 
definitions (|3S|) — (pp|) of the detector ideality can still be 
used. 

Let us emphasize that for the Bayesian description of 
the continuous measurement of a single qubit, the detec- 
tor in our model is completely described by six quantities: 
dc output current ("operating point") Iq, response A J, 
output noise Ss, ensemble dephasing rate Ts, correla- 
tion magnitude K, and the induced qubit energy shift Ae 
(the shift Ae can also have a classical contribution due 
to shift of the detector operating point). The quantities 
Sji, Ts, and K are analogous to the output, input, and 
cross-correlation noise terms, which are usually used for 
the description of a classical amplifier (and similarly foe 
the description of a quantum amplifier - see, e.g. Ref.O 
and references therein). The induced energy shift Ae is 
somewhat similar to the effect of a finite amplifier input 
impedance onto the previous stage parameters. 



VI. DETECTOR MEASURING ENTANGLED 
QUBITS 

Finally, let us consider the case when the detector is 
coupled to N arbitrarily entangled and arbitrary inter- 



acting qubits. Following Ref.0 we introduce the mea- 
surement basis consisting of 2 N states \i) and up to 2 N 
different levels Ii of the detector average current (some 
average currents can coincide). For an ideal symmet- 
ric detector and "frozen" qubits, H q b s = 0, where H q b s 
is the Hamiltonian of intrinsic evolution of the qubits, 
the evolution of the density matrix p of qubits due to 
measurement is described by simple "Quantum Bayes" 
equatior£j'c2l 



Pij(r) = Pij{0) 



y/Pj(r)Pj(T) 

EkPkk(o)p k (ry 



(51) 



where Pi(r) is the probability of obtaining particular 
measurement result (Id in this case) assuming state \i) of 
the measured system: 



Pi(r) 



:exp[ = _ (52) 



(we again assume that the currents Ii do not differ much, 
and therefore the detector noise £0 is state- independent). 

To take into account additional output noise £i(t) with 
spectral density Si , we have to perform the averaging of 
the density matrix pij over x = I — Id with the weight 
factor 



(53) 



p(x) = P{x)/ / P(x)dx, 



( x 2 \ r (i-x-lk) 2 



This procedure (without account of other noise sources) 
will lead to the equation presented in Ref.Cj and corre- 
sponds to the detector ideality r\ = So /(So + Si), similar 
to the one-qubit case. 

For the classical backaction noise which is 100% cor- 
related with £i(i), we should take into account that 
it can be coupled differently to different qubits. Let 
us assume that the energy of each state |z) is affected 
by the classical backaction noise, proportional to £1, so 
that Si — > Si + Aj£i (t), where Ai are arbitrary con- 
stants. Then Eq. (^lj) should be multiplied by the factor 
exp[i(A, — Ai)xr /K\ (in the one-qubit case the previously 
defined A would correspond to A2 — A\). 

Similarly, to take into account the possible asymmetry 
of the quantum backaction noise, let us assume that each 
electron tunneling through the ideal part of the detec- 
tor shifts the phases corresponding to states |i) (differ- 
ently for different states), that leads to the extra factor 
exp[i(0j - 9j)T d T/h] in Eq. @. 

The uncorrelated classical noise £3 (t) is also assumed to 
be coupled differently to the states \i), so that Si — > Si + 
<7j£3 (t), where gi are some constants. The averaging over 
noise £3 is simple and leads to the extra factor exp[— (gi — 
9j ) 2 S 3 T/4h 2 } in Eq. (||). 

Taking into account the effect of 9i , averaging over the 
noise £i(t) (and fully correlated backaction noise) and 
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£a(t), differentiating equations over time, and adding the 
intrinsic evolution of qubits, we finally obtain the follow- 
ing equation in the Stratonovich form: 



i Ae ■ 

Pij = -r[H q bs, p]ij + i—^-Pij + iKij[I(t) 



Ii + h 



\pij 



3 : s- ^2 



Pkk 



h + h 



}(ij-h) 



Jij Pij i 



(55) 



where the first term describes the intrinsic evolution of 
qubits due to H q t, s , Ae^ = — ®j)(Ii + Ij)/% is the 
effective energy shift due to detector asymmetry, Ss = 
So + Si is the total output noise, 



(Aj - A t )S 1 + (0i - Bj)S 
hSv 



(56) 



is the correlation factor between output and backaction 
noises, and the dephasing rate is 



(Ij-IjfSx ( 9i - 9j ) 2 S 3 
4S Sj; Afi 2 



(57) 



Notice that there are obviously no dephasing terms for 
diagonal matrix elements. Also notice that Eq. ( |55| ) is ap- 
plicable to both linear and ponlinear detectors, including 
purely quadratic detectors,cil as long as the condition of 
weak response is satisfied. 

Translating Eq. (p5l) into Ito form, we obtain 



Pij - -j\tiqbs,P\i 3 + 1— ^ 



-Pij 



\K ij [I{t)-Y,PkkI l 



-pij-q-\i{t) - y^,pkkh](ii + ij - zy^^Pkkh) 

s k k 

~^ijPiji 



k\ Pij 



(58) 



where the ensemble dephasing Tij is related to the single 
system dephasing 7^ as 



r = 

13 



{h-Ijf K*S* 
45 s 4 

and in our particular case is equal to 
^ _ {h-Ijf , (Aj-AtfSx 



' lij 



(59) 



0j) 2 S O 



45 

(ffj - gj) 2 5 3 

An 2 



Ah 2 



An 2 



(60) 



Since the combination I(t) — ^ k Pkklk m Eq. (58) is a 
pure noise because of the relation 



*(t) = X)p**-r* + &(*) + &(*)» 



(61) 



the ensemble averaged evolution is described by the re- 
duced equation 



i Ae.- 

Pij = --^[Hqbs,p]ij + '^—^-Pij - 1 ,.,/'-,• 



(62) 



Because of the reciprocity, it is natural to assume that 
the backaction couplings Aj — Ai, 9i — 8i, and gi — gj 
are proportional to the signal coupling Ii — Ij, so that 
Aj - Ai= (Ii - I )a, 6i - 9j = - Ij)Q, and g l - gj = 
(Ii — Ij)g. (Actually, this assumption implies detector 
linearity and also that all interactions with qubits occur 
via one "port of entry" . It is not valid, for example, 
when several geometrical parts of the detector interact 
with qubits in different ways.) With this assumption 
the parameters Ae^ 
equations ([[HI) and (|58j) become 



Kij, jij, and used in evolution 



Asi 



Kij = 



= e(i*-i?)/2, 

aS\ + OSq 



TiS-e 



^ = (h-I,) 2 [ 



Si 



Ii), 



ASqSt 



(a - Q^SqSx 
4^ 2 5 E 



' 73,' 



■ ij - (Ii Ij) ( AS(j 



1 a 2 ^ e 2 s* 



Ah 2 



Ah 2 



73,»), 



(63) 
(64) 

(65) 

(66) 



where 73 jT! = gS^/ATi 2 . Notice that there will be no 
dephasing between states \i) and |j) if the detector is 
equally coupled to these states, ij — Ij. 

The detector ideality in this case can be characterized 
by a single number (or few numbers for different defini- 
tions), which does not depend on the state of the mea- 
sured system. Extending the definitions (^J), fl39|), and 
( |40| ) discussed in previous Sections, the detector ideality 
can be characterized by the parameter combinations 



1/45 E _ 1/4S S + K 2 n S^/A 
V = -F > V = - 



r E , r 



r E> , 



V2 = 



1/45! 



i\ 



(67) 



where K n = (aSx + OS )/hSz and r S) „ = l/45 



a 2 S x lAh 2 + Q 2 S /Ah 



■ 73,-. 



In the case a = 9 = all 



definitions of ideality coincide and the evolution equation 
( p5| ) reduces to the equation derived in Ref.c9. In the case 
of finite a and/or 0, more natural definitions are fj and 
772 (again, fj 2 is the total energy sensitivity in units of 
H/2). However, ideality rj can also be a useful parameter, 
for example, if there is no way to control the degree of 
freedom affected by the backaction noise 6£o + a£i, and 
therefore the corresponding dephasing cannot be reduced 
by a feedback procedure. 



VII. CONCLUSION 

In this paper we have analyzed the process of continu- 
ous measurement of a solid state qubit by a nonideal solid 
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state detector. We have considered the phenomenologi- 
cal model of the detector (Fig. [I]) consisting of an ideal 
(quantum-limited) part and classical noise sources which 
contribute to the output (£1) and backaction (£2 + £3) 
noises. The possible correlation between classical out- 
put and backaction noise sources is taken into account 
by separating the backaction noise into a contribution 
£2(4) Mly correlated with output noise and the mi- 
correlated contribution £3 (t) . For the description of the 
ideal PEWjnW have started with the Bayesian equations 
of Refs.uu'EJ and then used the model of an asymmet- 
rically coupled ideal detector developed in Ref.E3. The 
asymmetric coupling changes the self-consistent energy 
difference between two qubit states. Also, this change 
fluctuates in time and the fluctuations are correlated with 
the output noise, thus producing an effect similar to the 
correlation of classical noises. 

The main result of the paper for the one-qubit case is 
the derivation of evolution equations (|30|)-(|3l|) and (|34|)- 
( |35| ) in Stratonovich and Ito form, respectively (the qubit 
energy parameter e should be treated self-consistently, as 
discussed in section V). In these equations the detector 
is characterized by the total output noise Ss, induced 
ensemble qubit decoherence rate Ts, and the total cor- 
relation K [see Eq. (|m])] between output and backaction 
noises, so that the detector separation into the quantum 
part and extra noises is irrelevant. (Notice that these 
three quantities arc the counterparts of output, input, 
and cross-correlation noise terms used for the descrip- 
tion of a classical amplifier; the induced qubit energy 
shift Ae is somewhat analogous to a backaction due to 
finite input impedance.) The relation between ensem- 
ble and single qubit decoherence rates is given by Eq. 
(|32|), which leads to the fundamental limitation ( |38|) for 
the ensemble decoherence rate. The discussed definitions 
of the detector ideality [see Eqs. (p^), (39), and (f4(i|)] are 
various combinations of the single qubit decoherence rate, 
ensemble decoherence, and the "information acquisition" 
rate (A/) 2 /45s. A 100% ideal detector corresponds to 
the absence of single qubit decoherence. 

The theory developed for a single qubit measurement is 
generalized to the case of entangled qubits in Section VI. 
The evolution equation is given by Eqs. ( |55| ) and (58), 
while the relation between ensemble and single system 
decoherence rates is given by Eq. ([59]) . 

The author thanks R. Ruskov for permission to use his 
result on ideal asymmetric detector prior to publication 
and for critical reading of the manuscript. The work 
was supported by NSA and ARDA under ARO grant 
DAAD19-01-1-0491. 
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